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Superfluidity and Bose-Einstein condensation in a dipolar Bose gas with weak disorder 


Abdelaali Boudjemaa 

Department of Physics, Faculty of Sciences, Hassiba Benbouali 
University of Chief P.O. Box 151, 02000, Ouled Fares, Chief, Algeria. 

We investigate the properties of a three-dimensional homogeneous dipolar Bose gas in a weak 
random potential with a Gaussian correlation function at finite temperature. Using the Bogoliubov 
theory (beyond the mean field), we calculate the superfluid and the condensate fractions in terms of 
the interaction strength on the one hand and in terms of the width and the strength of the disorder 
on the other. The influence of the disordered potential on the second order correlation function, the 
ground state energy and the chemical potential is also analyzed. We find that for fixed strength and 
correlation length of the disorder potential, the dipole-dipole interaction leads to modify both the 
condensate and the superfluid fractions. We show that for a strong disorder strength the condensed 
fraction becomes larger than the superfluid fraction. We discuss the effect of the trapping potential 
on a disordered dipolar Bose in the regime of large number of particles. 
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I. INTRODUCTION 


Recently, the study of Bose-Einstein condensate (BEC) 
in disordered potentials, known as ’’dirty boson problem” 
[H, has paid a special attention. Experimentally, it was 
first studied in the context of the motion of superfluid 
helium in porous Vycor glass 0. The problem of bo¬ 
son localization, the superfluid-insulator transition, and 
the nature of elementary excitations in lattice models 
with random site potentials characterized by Hubbard or 
equivalent models have been the object of several theo¬ 
retical investigations [1, Q and Monte Carlo numerical 
simulations [a-Q- Such random potentials suppress or 
may even completely destroy the long-range order re¬ 
lated to BEC Q leading to the transition to the so- 
called Bose glass phase [l|. This state is characterized 
by a finite compressibility, the absence of a gap in the 
single particle spectrum, and a nonvanishing density of 
states at zero energy. Bose systems with disorder have 
also been addressed in the uniform system without lat¬ 
tices i.e. continuum regime using the Bogoliubov model 
im, hydrodynamic approximation [TTI| . Beliaev-Popov 
diagrammatic technique [l^ [iJ . mean-field theory [TJ , 
and Monte Carlo simulations [^. The main finding of 
the above studies is that both BEC and superfluidity 
are depressed due the competition between the two-body 
repulsive interaction and the random potential. The ran¬ 
dom potential can be created using different techniques 
for example, laser speckles which are produced by focus¬ 
ing the laser beam on a glass plate (In . , wire traps 
represent magnetic traps on atomic chips [l9l| . Also, deep 
[20l| and incommensurable lattices provide a useful ran¬ 
dom medium [T^. 

On the other hand, there has been immense progress in 
creating BECs with strong magnetic dipoles [2ll-l^ , and 
heteronuclear molecules with strong electric dipoles ^dt - 
(e.g.. RbCs) as well as Rydberg atoms in electric fields 
. What renders such systems particularly fascinating 
is that the atoms interact via a dipole-dipole interaction 
(DDI) that is both long-ranged and anisotropic which 


leads to the observation of novel phenomena in ultracold 
atomic gases [UIH, [s^ . Recently, uniform dipolar Bose 
gas with a Gaussian disorder correlation function [sij 
and a Lorentzian (3^ , have been investigated using mean 
field theory at zero temperature. Most recently, a delta- 
correlated disorder is examined by the same group [s^ 
employing the Bogoliubov theory at finite temperatures. 
The main result emerging from these explorations is that 
the anisotropy of the two-particle DDI is passed on to 
both BEC and the superfluid fraction. 

In this paper, we study the properties of a three- 
dimensional (3D) homogeneous dipolar Bose gas in a 
weak rairdom potential with Gaussian correlation by 
employing a perturbation theory at fiirite temperature. 
The Gaussian-correlated disorder is defined by both the 
strength and the correlation length of the disorder (see 
Ea. (ll5|) ) which permit us to well control the interplay be¬ 
tween the disorder correlation length and the interaction 
strength. Therefore, this makes our model completely 
different from that of Ref.fs^ and thus, our results are 
also different in particular at zero temperature. Among 
the advantages of the Gaussian-correlated disorder po¬ 
tential is that its presence in ID Bose gas with short- 
range potential can undergo a finite-temperature phase 
transition between two distinct states: fluid and insulator 
[33| . In 3D case, BEC with a weak, Gaussian-correlated 
disordered potential can constitute a good model to de¬ 
scribe the Bose glass phase M- Moreover, we investi¬ 
gate, for the first time, the behavior of the anomalous 
density of a dirty BEC with DDI. This quantity which 
grows with interactions and vanishes in noninteracting 
systems [ssl - ls^ . is important to fully understand the in¬ 
terplay of disorder and interactions. The anomalous den¬ 
sity is often neglected in the literature without any physi¬ 
cal or mathematical justificatioirs. Indeed, it is usually of 
order or even larger than the noncondensed density (see 
Eg. (ITU) ) regardless of the system in which it is studied 
namely BEC with pure contact interaction Jp®, BEC- 
impurity mixtures (^. 1^ or dipolar BEC |43l| . Moreover, 
the absence of the anomalous averages renders the sys- 
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tern unstable. We show, in addition, how the anisotropy 
of the DDI enhances quantum, thermal and disorder fluc¬ 
tuations as well as the superfluid fraction. We find that 
the thermal fluctuations’ contribution to the condensate 
depletion and the superfluid fraction coincide with those 
obtained in [s^ for a delta-correlated disorder. This is in 
fact natural since we assume that disorder, quantum and 
thermal fluctuations are too small that the Hamiltonian 
of the system can be expanded in both leading order of 
disorder plus thermal and quantum fluctuations. There¬ 
fore, disorder and thermal fluctuations are additive and, 
thus, independent from each other. 

The rest of the paper is organized as follows. In sec¬ 
tion im we describe our model of the dipolar dilute Bose 
gas in a general random potential. In section mi we 
derive analytical expressions for the condensate fluctu¬ 
ations, the second order correlation function and some 
thermodynamic quantities for external random poten¬ 
tial with Gaussian correlation at finite temperature and 
compare our finding with the mean field results of Ref 
We show that the competition between both con¬ 
tact interaction-disorder and DDI-disorder leads to en¬ 
hance the condensate depletion, the anomalous density, 
disorder fluctuation and the ground state energy. We es¬ 
tablish the validity criteria of the Bogoliubov theory of 
disordered dipolar systems. In section El the superfluid 
fraction is obtained and its characteristics are discussed. 
In section El we extend our results to the case of a har¬ 
monically trapped gas by making use of the local density 
approximation. Finally, our conclusions remain in sec- 

tioidl 


II. THE MODEL 

We consider the effects of an external random field 
on a dilute 3D dipolar Bose gas with dipoles oriented 
perpendicularly to the plane. The Hamiltonian of the 
system is written as: 

H = y d^r^’^(r) ^-^-A -|- U{r)^ ■!/>(r) 

-I-i / d^r f d^r'ip''‘{r')V{r — (1) 


where the coupling constant Cdd is for particles 

having a permanent magnetic dipole moment M {Mq is 
the magnetic permeability in vacuum) and d^ /cq for par¬ 
ticles having a permanent electric dipole d (eo is the per¬ 
mittivity of vacuum), m is the particle mass, and 0 is 
the angle between the relative position of the particles f 
and the direction of the dipole. The characteristic dipole- 
dipole distance can be defined as r* = mCdd/^T^fi^ ■ The 
disorder potential is described by vanishing ensemble av¬ 
erages {U{r)) = 0 and a finite correlation of the form 
{U{r)U{r')) = R(r,r'). 

Passing to the Fourier transform and working in the 
momentum space, the Hamiltonian m takes the form: 

TT \ ^ ^ .V t ^ X ^TT .. t 

(3) 

k k,p 

^ ^ ^ (p)®k4q®k-q®k+p°'k-p, 

k,q,p 

where B is a quantization volume, and the interaction 
potential in momentum space is given by 

/(k) = g[l-h edd(3cos^6»fe - 1)], (4) 

here Cdd = Cdd I^9 is the dimensionless relative strength 
which describes the interplay between the DDI and short- 
range interactions. 

Assuming the weakly interacting regime where r* <C ^ 
with ^ = h/y/mgn being the healing length and n is the 
total density, we may use the Bogoliubov approach. Ap¬ 
plying the inhomogeneous Bogoliubov transformations 

i 

dk = Mfedk-'yfedL^-/3k, dl = Mfe6k-'yfcS_k-/3k, (5) 

where and Sk are operators of elementary excita¬ 
tions. The Bogoliubov functi ons u^, Vk a re expr essed in 
a standard way: Uk,Vk = {\/ek/Ek ± with 

Ek = h^k'^/2m is the energy of free particle, and 

P^ = J^^Uk. (6) 

The Bogoliubov excitations energy is given by 


where and ip denote, respectively the usual creation 
and annihilation field operators, the interaction potential 
V{y — r') = gSpe — r') -|- — rOi 9 = 47r?i^a/m corre¬ 

sponds to the short-range part of the interaction and is 
parametrized by the scattering length as a. It is impor¬ 
tant to mention here that the contact potential exhibits 
ultraviolet divergences in both 2D and 3D geometries. 
Common ways to cure these inconsistencies is the use of 
either renormalization or a pseudopotential. 

On the other hand, the dipole-dipole component reads 

.. / X Cdd 1 - 3cos^6> 


Ek = ^El + 2tiod{0)Ek, (7) 

where Hodid) = n lim /(k) is the zeroth order chemical 
fc—>0 

potential. 

Importantly, the spectrum m is independent of the ran¬ 
dom potential. This independence holds in fact only in 
zeroth order in perturbation theory; conversely, higher 
order calculations render the spectrum dependent on 
the random potential due to the contribution of the 
anomalous terms (see below). For /c —>■ 0, the excita¬ 
tions are sound waves Ek = hcsd{d)k, where Csd{d) = 
Cs^/l + edd(3cos^ 9—1) with Cs = \fgnjrn is the sound 


47r 
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velocity without DDL Due to the anisotropy of the dipo¬ 
lar interaction, the sound velocity acquires a dependence 
on the propagation direction, which is fixed by the angle 
9 between the propagation direction and the dipolar ori¬ 
entation. This angular dependence of the sound velocity 
has been confirmed experimentally (^ . 

Therefore, the diagonal form of the Hamiltonian of the 
dirty dipolar Bose gas can be written as 

i? = + (8) 

k 


where E = Eqj, + SE + Ej^, 

Eodid) = 9-od{.0)N/2 with N being the total number of 
particles. 

5E=\Y.\-^k-Eu-nf{\^)\, (9) 

k 

is the ground-state energy correction due to qunatum 
fluctuations. 

Er =-^n{\Uk\^)^ =-^nRk^, (10) 

k k 


gives the correction to the ground-state energy due to the 
external random potential. 

The noncondensed and the anomalous densities are de¬ 
fined as h = X)k(®k®k) and to = ^j,(dkd-k), respec¬ 
tively. Then invoking for the operators Ok the trans¬ 
formation ®, setting (fo^fok) = i^k'k-^k and putting 
the rest of the expectation values equal to zero, where 
Nk = [exp(£k/T) — 1]“^ are occupation numbers for 
the excitations. As we work in the thermodynamic 
limit, the sum over k can be replaced by the integral 
= V f d?k/[2 ttY and using the fact that 2A^(a;)-|-l = 
coth(a:/2), we obtain: 



cPk Ek + /(k)n 
(27r)3 Ek 




+nR, 


( 11 ) 


and 


TO = — 


1 

2 


cPk fCk)n 
(27r)3 sk 


coth 



+ nR. 


( 12 ) 


The contribution of the random potential comes through 
the last terms in Eqs (HU and (HU- These terms are 
defined as 


riR = ^^{\l3k?) = n f 

k 


(2^)3 '=4' 


(13) 


Expressions (HU) and HU must satisfy the equality 


flkiflk + 1) - 


1 

dsinh^ (efe/2T) 


+ nR 



Equation (|14l) clearly shows that rh is larger than h at low 
temperature irrespective of the presence of an external 
random potential or not. So the omission of the anoma¬ 
lous density in this situation is principally an unjustified 
approximation and wrong from the mathematical point 
of view [iU SO, Si • 


III. FLUCTUATIONS AND THERMODYNAMIC 
QUANTITIES 


To proceed further in practical calculations, we must 
specify the type of random potential. Eor this purpose, 
we take the case of a spatially decaying disorder cor¬ 
relation R{r). Therefore, in what follows, we restrict 
ourselves to the case of a Gaussian correlation with the 
Eourier transform S, IsTjl 

i?(k) = (15) 


where R with dimension (energy) ^ x (length)3 and a 
characterize the strength and the correlation length of the 
disorder, respectively. Equation (ITSl) makes the macro¬ 
scopic wave function of BEC not sensitive to disorder in 
and between pores, but instead depends on the disorder 
averaged over the coherence length. Hence the ensemble- 
averaged system can become nearly uniform M- 

Substituting (fT^ in Eq. (|T^ , we obtain for the conden¬ 
sate fluctuation due to the disordered potential 


nR 


tn^R 

87r3/2?i^ 


a 


(16) 


where 


, . d9 sin 9E(a) 

n{edd,ct)= / - , - , 

Jo 2y^l -I- edd(3cos2 9k - 1) 


(17) 


is depicted in EigHJ 

The function E{a) = e^°‘{4,a + 1) [l — erf(-\/2a)] — 
2\j2aj'K with a = cr^[l -I- edd(3cos^d — 1)]/^^, has the 
following asymptotics for small a: E{a) = 1 —4\/2a/7r-|- 
6a — {2)2/‘3)ay/2apK + lOa^ -I- O Equation. (IT6l) 

is in good a gree ment with that obtained using the mean 
field theory im. 

We observe from EiglUthat for a small disorder corre¬ 
lation length (T < ^, the contribution of the DDI on the 
disorder fluctuation is not important whereas in the case 
of cr > ^, the DDIs tend to enhance the disorder fluctua¬ 
tion. 

On the other hand, for a/^ —>■ 0, we get from Eg. dTTl) 
that h{edd,6) = Q-ii^dd)- Thus, the disorder fluc¬ 
tuation (1161) becomes identical to that obtained in 3D 
dipolar BEC with delta-correlated disorder [sj nR = 
{m‘^R/8Tr^^'^h'^)y/n/aQ-i{edd), where the contribution 
of the DDI is expressed by the functions Qj{edd) = 

(1 - eddy^^ 2 Fi (-§, I; , where 2 E 1 is the hyper¬ 

geometric function. Note that functions Qj{edd) attain 
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FIG. 1. (Color online) Behavior of the disorder function 
h{edd, O'/^) from Ea. (ll7ll . as a function of ct/^. Black line; 
Cdd = 0 (pure contact interaction), blue line: edd = 0.15 (Cr 
atoms), red line: edd = 0.38 (Er atoms), green line: edd = 0.7 
and brown line: edd = 1. 


their maximal values for edd ~ 1 and become imaginary 
for e..>llilEi. 

For cr/^ —>• 0 and edd = 0, we read off from Ea. fTTll 
that one obtains h{edd,Oi) —>■ 1. Therefore, we should re¬ 
produce the Huang and Meng result @ for the disorder 
fluctuation in this limit. 

Upon calculating the integral in Eq. dm, we get for the 
noncondensate depletion 


n 

n 


8 na^ 


“o \/) Q-li^dd) 




+ 2ttR'\I —^h{edd, a), 


(18) 


where R' = R/g^n is a dimensionless disorder strength. 
The condensed fraction can be calculated employing 
Tic/n = 1 — fi/n. 

The integral in Eo. lfT^ is ultraviolet divergent. This 
divergence is well-known to be unphysical, since it is 
caused by the usage of the contact interaction potential 
as we have mentioned in SecUll A general way of treating 
such integrals is as follows. First, one restricts to asymp¬ 
totically weak coupling and introduces the Beliaev-type 
second order coupling constant [i^ 


the anomalous fraction turns out to be given as 


m 

n 


na^ ^ , 

=^]J — Qsi^dd) 


2 Ina^ 

3 V ~ 



2 

Q—1 (cdd) 


+ 2nR\l^h{edd,a). 


( 20 ) 


The leading terms in Eas. lfTSl) and (1^ represent the 
qunatum fluctuationji^. The subleading terms which 
represent the thermal fluctuationji^, are calculated at 
temperatures T <C gn, where the main contribution to 
integrals (HU and m comes from the region of small 
momenta (e^ = fiCsdk). The situation is quite different 
at higher temperatures i.e. T ^ gn^ where the main con¬ 
tribution to integrals (ED and ED comes from the single 
particle excitations. Hence, the thermal contribution of h 
becomes identical to the density of noncondensed atoms 
in an ideal Bose gas [i^, while the thermal contribution 
of m tends to zero since the gas is completely thermal- 
ized in this range of temperature [H, 113, . The last 

terms in ED and (HOD describe the effect of disorder on 
the noncondensed and on the anomalous densities. 

Equation (EOD clearly shows that at zero temperature, 
the anomalous density is three times larger than the non¬ 
condensed density for any range of the dipolar interaction 
as well as for any value of the strength and the correlation 
length of the disorder as it has been anticipated above. 
Moreover, rh changes its sign with increasing tempera¬ 
ture in agreement with uniform Bose gas with a pure con¬ 
tact interaction [i^. For edd = 0, Qj{edd) = 1 and thus, 
Eqs. ED and (HOD reproduce the short-range interaction 
results. Furthermore, the DDI enhances the condensate 
depletion and the anomalous fraction for increasing edd- 

The condensate fluctuations manifest into the second- 
order correlation function as [i^ 


G^^\r) = (i^’^(r)V)’^(r)V)(r)'0(r)) 

= nl + + 2h^ -I- dhrzc -b 2mnc. (21) 

Equation m is obtained using Wicks theorem. Inserting 
then Eas. lfTBl) . (fTSl) and (001) into ED, we obtain 

. 64 jna^ ^ \ ^ 

— 2 “ = 1 + tV- Qai^dd) + 8 ttR \ - h{edd,a). 

n 6 \ IT V TT 

( 22 ) 

This equation is accurate to the first order in h/ric and 
mine and shows how the correlation function depends on 
the interaction parameter edd and the disorder length a. 

The energy shift due to the interaction and the quan¬ 
tum fluctuations ([9]) are ultraviolet divergent. The dif¬ 
ficulty is overcome if one takes into account the second- 
order correction to the coupling constant ED- A 
straightforward calculation yields]^ 


Mk) =/(k) - (19) (23) 

However, the energy shift ED due to the external ran- 
After the subtraction of the ultraviolet divergent part, dom potential ()15p is not divergent and it can be evalu- 
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ated as 


Er 


= WttR' 


nw- 


{^dd: Ck), 


where Eq = Ngn/2, and 


(24) 


X _ 

hi{€dd,a) = - d6'smv^l + edd(3cos2 Ok - l)Ei{a), 

(25) 


is displayed in FiglU 

The function Ei{a) = e^“[l — erf(-\/2a)] — y^l/27ra has 
the asymptotics behavior for small a\ Ei{a.) = 1 — (4 + 
'n)\/aj^-K + 2 a — ( 8 / 3 )-\/ 2 / 7 ra^/^ + 2 a^ + O (a®/^). 

As is seen from FigI5] that for a < 2^, the energy cor¬ 
rection due to the disorder effect (1551) is negative which 
leads to lower the total energy of the system. Note that 
this result still valid for any value of tdd < 1- Another 
important remark is that the energy decreases with in¬ 
creasing edd- 

For a condensate with a pure contact interaction 
{Q 5 {^dd = 0 ) = 1 ) and in the absence of disordered po¬ 
tential {R = 0), the obtained energy excellently agrees 
with the seminal Lee-Huang-Yang quantum corrected 
ground state energy (4^ . 

For, a/^ —>■ 0, the energy shift due to the external ran¬ 
dom potential m becomes ultraviolet divergent. Again, 
by introducing the renormalized coupling constant dUD 
one gets: Er/Eq = IdirR'y/na^/nQi^Cdd) which well 
coincides with the result obtained with delta-correlated 
disorder of Ref @ . 



FIG. 2. Behavior of the disorder energy function hi{edd, c/^), 
from Eq. dlSl) as a function of a/^ for same values of edd as 
in Figd] 

The chemical potential can be obtained easily through 

H = dE/dN or via /r = /rg + 5^, _ where dfi = 

J2kfi^){i’k{vk-Uk)] = Z)k/(k)(^ + w) One 


can also calculate the shift of the sound velocity utilizing 
mc^ = ndfi/dn [i^ . 

The Bogoliubov approach assumes that fluctuations 
should be small. We thus conclude from Eqs. m 
and (l20l) that at T = 0, the validity of the Bogoli¬ 
ubov theory requires inequalities Vna^Q^^edd) ^ 1 and 
R'Vna^h{edd, a) -C 1. For R' = 0, this parameter differs 
only by the factor Q 3 {edd) from the universal small pa¬ 
rameter of the theory, VruE <C 1, in the absence of DDL 
At T ^ gn, t he Bogoliubov theory requires the condition 
(Tlgn)Vna^Q-i{edd) ^ 1- The appearance of the extra 
factor (T /gn) originates from the thermal fluctuations 
corrections. 


IV. SUPERFLUID FRACTION 

The superfluid fraction Us/n can be found from the 
normal fraction n„/n which is determined by the trans¬ 
verse current-current correlator Us/n = 1 — Un/n. We 
apply a Galilean boost with the total momentum of the 
moving system P = mvinVs + ^nVn), where Vg denotes 
the superfluid velocity and and Vn = u —Vg is the normal 
fluid velocity with u being a boost velocity [s^- Using 
da, the superfluid fraction is then written 

f (Pk nRkkikj 

n J (27r)3 2m - 2n/(k)]2 


2 

f (Pk 

fP kikj 

Tnj 

(27r)3 

.2m4sinh^(£fc/2r)_ 


It is worth noticing that if in expression (1261) rh were 
omitted, then the superfluid fraction would be divergent 
signaling the importance of the anomalous density on the 
occurrence of the superfluidity in Bose gases [H, . 

Equation. (ES) yields a superfluid density that depends 
on the direction of the superfluid motion with respect to 
the orientation of the dipoles. In the parallel direction 
and at low temperatures where £k = TiCsdk, the superfluid 
fraction reads 


= 1 - 47rR' 


na 

TT 


3 p_.2^4 

/i"(edd,a) - 77- ^^^l^iedd), 


where the function 


h''''{edd,a) = [ 
Jo 


dte- 


A5mnfPc^ 


sin 6 cos^ OR {a) 


0 2y^l -I- edd(3cos2 9-1) 


(27) 

(28) 


is decreasing with increasing edd and vanishing for large 
tj/^ as is depicted in FiglSla. And the functions 

Q^ji^dd) = 7(1 - ^ddy^'^iFi (-§, 7 ^^), have the 

properties Q^j{edd = 0) = 1/3 and become imaginary for 
Edd > 1 (see Figll]). Therefore, Eg. (1571) reveals that DDI 
effects are more significant for condensate fracion m 
than for the parallel superfluid fraction. 



















h^{edd,cr/^) h"(edd,(T/^) 
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FIG. 3. (Color online) Behavior of the disorder functions 
h^^{edd,cr/^) (a) and h^{edd,o/£) (b), as a function of cr/^ for 
same values of tdd as in Fig[T] 


Again at low temperatures, the perpendicular direction 
of the superfluid fraction (E51) takes the form 



n 


= 1 - 2ttR\/ ^^h-^{edd, a) 


tt^T^ 

Abmnh^cl 


Q-^i^dd), 

(29) 


where the functions 


^ (cdd; tr) 



de 


sin 0(1 — cos^ 0)T{a) 
2\/l + edd(3cos2 0-1) 


= h{edd,a) - h}\edd,a), 


(30) 


and Qf{edd) = Qj{^dd) - Q^kedd), are displayed in 
Figs. (Ob), and (O, respectively. 



FIG. 4. (Color online) Behavior of thermal functions q|| 
(Green line) and Qf (Orange line), as a function of edd- 

Expressions (EH) and (12911 constitute a natural exten¬ 
sion of those obtained in since they contain the tem¬ 
perature correction (third terms). At T ^ ng, it is ev¬ 
ident that both thermal terms of rig coincide with the 
noncondensed density of an ideal Bose gas. Figure (|4]) 
shows that the thermal contribution of is smaller than 
that of nl for edd < 0.5, while the situation is inverted 
for edd > 0.5. 

For a/^ —7> 0 and tdd = 0, both components of the super¬ 
fluid fraction (l?7l) and (1^ reduce to rig/n = 1 — du^/Sn, 
which well recove earlier results of Refs S El, m for 
isotropic contact interaction. For a/^ —>■ 0, we have 

h^^{edd,0) = Q-ii^dd) and h^{edd,0) = Q^i{edd)- Con¬ 
sequently, the disorder correction to superfluid fraction 
m becomes identical to that obtained in 3D dipo¬ 
lar BEC with delta-correlated disorder [s^. We should 
stress also that for increasing edd, h^Kedd,ct) decreases, 
whereas h{edd,Oi) increases for fixed a/^. Therefore, this 
reveals that there exists a critical value beyond which 
the system has the surprising property that the disorder- 
induced depletion of the parallel superfluid density is 
smaller than the condensate depletion even at T = 0. 
This can be attributed to the fact that the localized par¬ 
ticles cannot contribute to superfluidity and, hence, form 
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obstacles for the superfluid flow. For a large disorder 
correlation length i.e. a ^ decreases indicating 

that the locally condensed particles are localized in the 
respective minima for the disorder potential for a finite 
localization time [s^- 

The superfluid fraction can be either larger or smaller 
than the condensate fraction Uc/n, depending on temper¬ 
ature, the strength of interactions, and on the strength 
of disorder. Increasing R' leads to the simultaneous dis¬ 
appearance of the superfluid and condensate fractions. 
So for any value of na^ and of Cdd there exists a critical 
strength of disorder 

d || ^ ^ Qsi^dd} 7*^11 

7rh»(€dd,a)’ “ 7rh^(€dd,a)’ 

for which n\/n < nc/n and jn < ric/n. 

When cr/^ -j. 0, i?!! = (4/7r)[Q3(edd)/Qii(edd)] and 
= {2/T:)[Q^(tdd)/Q-i{e-dd)\ ■ In the case of Er atoms 
(cdd = 0.38), rI ~ 6.74 and R^ ~ 0.78. For Cr atoms 
(cdd = 0.15), rI ~ 4.96 and R^ ~ 0.86. This clearly 
shows that i?!! is decreasing with Cdd, while R^ is in¬ 
creasing with Edd- 

Therefore, the Bogoliubov approach should satisfy the 
condition:!?' < Rc. However, it is not clear whether these 
results are still valid for R' > Rc in a range of densities 
where the difference between ns/n and Uc/n can be sig¬ 
nificant and hence, the system yields a transition to a 
new quantum regime. The response to these questions 
requires either a non-perturbative scheme or numerical 
Quantum Monte Carlo simulation, which are out of the 
scope of the present work. 


V. TRAPPED DIPOLAR BEG WITH WEAK 
DISORDER 

In this section, we discuss the case of a harmonically 
trapped dipolar Bose gas with weak disorder. Let us 
start by writing the hydrodynamic equations. These 
equations can be derived by factorizing the condensate 
wave function according to the Madelung transformation 
(/)(r,t) = Y^n(r, t) exp(—i5'(r, t)), where S is the phase of 
the order parameter, it is a real quantity related to the 
superfluid velocity by u = Ti/iriS/S. The continuity and 
Euler-like equations read, respectively: 

3ti 

^+V.(nu)=0, (32) 

dv \ A^Jn 1 9 

“a =+ r*” +" + 9 " + ®^- . 

(33) 

where Ayjnjyjn stands for the quantum pressure, ^dd = 
/ dv' Vd{v - r')|0(r')p and U{r) = Utraph) + U{r) with 
Utrap{i^) = rn{uj‘^x‘^ +ojyy‘^ +uj‘lz‘^)/2 being the trapping 
potential and ujj are the trapping frequencies. 


In a nonstationary situation, we consider small oscilla¬ 
tions (low density) for the density around its static solu¬ 
tion in the form n = no -|- (5n, where Sn/uQ <C 1. 

Shifting the phase by —fit/h, we then linearize Eg. (15^ 
with respect to Sn and VS' around the stationary solu¬ 
tion. The zero-order terms give for the chemical poten¬ 
tial: 

y. = -^^^^+U{r) + gno{r) + <^dd- (34) 

2m y/no 

Expanding the density and phase in the basis of the ex- 
citations S(r) = ^-f/2y'no(r)^ - h.c.] and 

i5n(r) = yno(r) {r)bk -f h.c.], one then finds that 

/^ = Mfe ± Vk obey the Bogoliubov-de Gennes equations 
(BdGE): 


-—A+U{r) +gno{r) +^dd- y 
Zm 


-—A+ll{r) + 3gno{r) + ^dd- fd 
2m 


/fc^W =efe/fe (r), 

(35) 

fk{'^)=^kfk{^)- 


Equations ([34|) and (IMll form a complete set to calcu¬ 
late the ground state (BEG) and excitations of the dirty 
dipolar Bose gas, from which one can compute all the 
properties of finite temperature, collective modes or time- 
dependent dipolar BEG with random potential. 

To derive analytical expressions for the physical quan¬ 
tities of interest such as the condensate depletion, the 
equation of state, the ground state energy and so on, one 
should use the local density approximation (LDA) (for 
more details on the applicability of this approximation 
on long-range interactions, see e.g. 0). The LDA or 
semiclassical approximation is applicable when the ex¬ 
ternal potential is sufficiently smooth, and requires that 
(r) are slowly varying functions of the position as well 
as U (r) is assumed to be a weak random potential. Em¬ 
ploying the LDA, the BdGE (l3^ take a simple algebraic 
form as in the homogeneous case. Therefore, the local 
Bogoliubov spectrum can be obtained in the usual way 
and reads £fe(r) = + 2no{r)fkEk. Then, we obtain 

for the condensate depletion, the anomalous fraction and 
the ground state energy the following relations, respec¬ 
tively: 


h(r) 

«o(i’) 


lm{v)a^ 


Qsi^dd) + 2ttR')J~- h{edd, cx) + ■ 

(36) 


«o(i’) 

and 


iW o,/«o(r)a3 


,, /no(r)a3 


Qsi^dd) + 2TrR'\l --- h{edd, a) + 

(37) 


TDi 'Ro{r)a^ 


Eo 


= UttR 


-hi{edd, a), 


(38) 




























where Eq = Ngno{r)/2. 

In the regime where the number of particles is very large 
(Thomas-Fermi limit), the density no remains extended 
(delocalized) for a weak U: no(r) = [g —U{t) — ^dd)/g- 
Equivalently, one can also use Eas. (l36l) and (l37)) to cal¬ 
culate the corrections to the equation of state and the 
local superfluid density. Importantly, we see that the 
disorder corrections remain anisotropic as in the homo¬ 
geneous case. 


VI. CONCLUSIONS 

In this paper, we have studied the properties of a dipo¬ 
lar Bose gas in the presence of a weak random poten¬ 
tial with a Gaussian correlation at finite temperature. 
Within the Bogoliubov approach, we have calculated the 
condensate fluctuation due to disorder, as well as the 
corresponding corrections to the noncondensed and the 
anomalous densities, second order correlation function, 
the ground state energy and the superfluid fraction in 
the homogeneous case. We have pointed out that the 
interplay of the DDI and the external random potential 
makes both the BEG and the superfluidity anisotropic. 
In particular, we have discussed the consequences of the 
superfluid density, which becomes a tensorial quantity as 
a linear response to the moving disorder. We have found 
that the presence of the DDI causes a weak anisotropy in 
the parallel component of the superfluid density while the 
perpendicular component becomes strongly anisotropic. 


We have also demonstrated that for a strong disor¬ 
der strength the system introduces an unusual quantum 
regime where the superfluid fraction is smaller than the 
condensate fraction. In addition, we have reproduced the 
expression of the condensate fluctuations and thermody¬ 
namic quantities obtained earlier within the treatment 
of Huang and Meng in the absence of the DDI and 


those obtained recently in Ref [33| in the limit cr/^ —>■ 0. 


Furthermore, we have discussed the validity criterion of 
the Bogoliubov approach in a disordered dipolar BEG at 
finite temperatures. By means of the LDA, we have gen¬ 
eralized our results to the case of a harmonically trapped 
gas. These results could be directly applied to check how 
quantum and disorder fluctuations can alter the time-of- 
flight expansion of a dirty-trapped dipolar Bose gas. 

Promising candidates for the experimental realization 
of such dirty dipolar BEGs are atomic species with highly 
magnetic dipolar interaction such as Dy (magnetic mo¬ 
ment lO^s) [ 2 ^ or polar heteronuclear molecules such as 
KRb (magnetic moment 0.6 Debye) [l^. 

Finally, an important step for future theoretical studies 
is to analyze the effects of an external random potential 
in a dilute 2D dipolar Bose gas near the roton minimum. 
In ID case, the investigation of a dirty dipolar BEG in an 
optical lattice is expected to exhibit even richer insights 
about the stability of the superglass state in Bose gases 

B- 
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